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Abstract 

We discuss the prospects of searching for the neutral Higgs bosons 
of the triplet model in central exclusive production at the LHC. A 
detailed Monte Carlo analysis is presented for six benchmark scenarios 
for the Higgs boson, these cover m^jo = 120, 150 GeV and doublet- 
triplet mixing of ch = 0.2, 0.5 or 0.8. We find that, for appropriate 
values of ch, an excellent Higgs mass measurement is possible for the 
neutral Higgs in the triplet model, and discuss how to distinguish the 
triplet model Higgs boson from the Higgs boson of the Standard Model. 

1 Introduction 

It is expected that the nature of electroweak symmetry breaking will be re- 
vealed by the LHC experiments in the near future. In the Standard Model 
(SM), the electroweak symmetry is spontaneously broken by a Higgs dou- 
blet, which contains a neutral scalar field that acquires a vacuum expec- 
tation value (VEV). However, several Higgs multiplets typically occur in 
extensions to the SM. In supersymmetric models, at least one additional 
doublet is required. In left-right symmetric models, triplets are added to 
naturally generate a small mass for the neutrinos. Although the new scalars 
do not always take part in the electroweak symmetry breaking, they affect 
the properties of the Higgs boson through mixing. 

Models with an extended Higgs sector typically contain charged scalars. 
A large number of studies [HE] have previously investigated the possibility of 
studying the doubly or singly charged components of higher representatioE0. 



^For example the discovery of the ZW''^ and/or W W H^^ vertices would 

serve as a direct proof of the non-standard structure of the Higgs sector (see e.g. [3]). 
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However the charged scalars may be considerably heavier than the hght 
neutral bosons. Therefore, it would be instructive to study the properties of 
the light neutral Higgs particles in order to reveal the manifestation of new 
representations [2]. 

Higgs triplets are an especially attractive possibility [5j. A tiny neutrino 
mass may indicate that the mass is being generated by the seesaw mechanism 
containing the coupling of neutrinos to the triplet. In addition, composite 
Higgs models contain several multiplets, including the triplet ones. Triplets 
also occur in the little Higgs models - see, for example, [6] and references 
therein. 

Determining that a new detected state is indeed a Higgs boson and dis- 
tinguishing it from the Higgs boson of the SM will be far from trivial. This 
task will require a comprehensive programme of precision Higgs measure- 
ments. In particular, it will be of utmost importance to determine the spin 
and CP properties of a new state and to measure precisely its mass, width 
and couplings. In this work, we suggest that the neutral sector of the triplet 
representation can be studied using the central exclusive production (CEP) 
mechanism (see, for example, [7]) if forward proton detectors are installed 
at ATLAS and/or CMS, (see [8]). 

The structure of the paper is as follows. In Sec. [21 we consider the general 
properties of models with higher representations. We then concentrate on 
the triplet representation in Sec. [3l we choose a benchmark model with the 
electroweak p-parameter equal to unity at tree-level, although the results 
are quite general. In Sec. HI we introduce the central exclusive production 
process. Finally, in Sec. El we present a detailed Monte Carlo analysis of the 
central exclusive production of a neutral Higgs boson in the triplet model 
for a selection of parameter choices. 



2 Models with general Higgs representations 

We start with the Standard Model gauge group SU(2)/,xU(l)y for the elec- 
troweak sector. The masses of the gauge bosons are then obtained from the 
kinetic part of Lagrangian, 

k i 

where are complex representations and are real ones. The covariant 
derivative is 

D^ = d^ + igW^T'=^ + ^g'B^, (2) 

where is the generator of SU(2) in the appropriate representation (with 
Tr(r'^r^) = ^(5"^) and Y is the U(l) hypercharge. Here W and B are the 
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SU(2) and U(l) gauge bosons respectively, and the mixing angle 9w of the 
Z boson and photon is obtained by diagonalizing the neutral sector. The 
W and Z boson masses are given by 

ml = {9' + 9")J:Ti^vl ml. = 9'T.T3^vl (3) 

i i 

where T^i is the isospin third component and Vi is the VEV of particle i. It 
is clear from Eq. ^ that the doublet VEV decreases when several repre- 
sentations obtain non-vanishing VEVs. Furthermore, since the left-handed 
fermions are in doublets, the charged fermions can only get their masses 
through the Higgs doublet representatiorH, ruf = y fVdoubiet , and the fermion 
Yukawa coupling, yf, must increase to produce the fermion masses. This, for 
example, leads to an enhancement in the production cross section for Higgs 
production via gluon fusion, where the dominant contribution is due to the 
top quark loop. A further enhancement is present in the branching ratio to 
fermion anti-fermion pairs. The possibility arises, therefore, of observing a 
very different prediction to that of the Standard Model. 

The higher Higgs representations are severely restricted by the elec- 
troweak p-parameter. The p-parameter in the Standard Model is defined 
by the ratio of the gauge boson masses, 

p = , (4) 

m| cos^ 6w ' 

which at tree level is exactly unity in the Standard Model. The radiative 
corrections to the «-parameter have been studied in the SM up to three- 
loop level [U] - [15)^1 . For rriH ~ 2m vi^, the correction to the p-parameter 
is p — 1 ~ — 0.00078+(4-loop and higher corrections). For heavier Higgs 
masses, the absolute value of the negative corrections increase. In a model 
with several scalar representations, whose neutral component develops a 
VEV, the p-parameter is given at tree level by [18j 



E^r, (T,(T, + 1)-T, 



3i)^i 



(5) 



Here Tj is the weak isospin and rj = 1/2(1) for real (complex) representa- 
tions (see Eq. ([7]) for examples). Finally, the p-parameter is experimentally 
constrained to be [19] . 

p-l = 0.0002 , (6) 



^For the neutral fermions this is not the case, since the Majorana masses can be 
generated through triplets. 

^An explicit formula for the one-loop correction Ap*-^' in p = Ap'^' + Ap''^' + Ap'^' is 
given in [l^, two- loop terms are given in [17], and three-loop corrections in [15j . 
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where the quoted errors are at 2a. As the loop-corrections to the /j-parameter 
in the SM are negative, it can be argued that the p-parameter favors either 
a light Higgs boson or models that result in positive corrections to p. 



3 Higgs bosons in a triplet model 

In order to fulfill the experimental constraint on the p-parameter in Eqn. Q, 
the triplet VEV has to be small. Using Eqs. (0) and © one finds that the 
upper limit for the triplet VEV is a few GeV's assuming that there is one 
triplet in addition to one doublet. An alternative method to satisfy the 
experimental constraint at tree-leveE| is to have representations which add 
up to p = 1. We discuss this next. 

3.1 A model with p = 1 

We consider the model developed by Georgi and Machacek [Ij, and further 
studied in |22| [23] in which additional representations are chosen in such 
a way that the tree-level value of p remains unity. The p-parameter is fixed 
to one by choosing one complex scalar doublet (i^y=i) and two triplets, one 
real (^y=o) and one complex (xy=2)- These can be written as 



The following sign conventions are chosen: (j) = —(0"^)*, x = (x^"^)*! 
X" = -(x+)*, i~ = -{i^Y, and = The VEVs of the neutral 

components of the Higgs fields are denoted by (x*^) = = b and {(j)^) = 
a/\/2- For doublet-triplet mixing, the standard notation is employed: 



*At one-loop level, one has to consider renormalization. It has been shown in i20^ that 
p ^ 1 at tree-level is acceptable in a real triplet model as the experimental measurement 
of p is satisfied after calculating higher order corrections. 






X -X, is 



V = Ai(TV</>V-4^;2)2 + A2(Trx^x-^4^')' 
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Figure 1: Masses of the Higgs bosons, mn, , divided by the unknown y/Xi 
couplings, as a function of ch- 



where rj/2 are the SU(2) generators in the doublet representation and ti in 
the triplet representation. 

As we are interested in this model mainly to illustrate the possibility of 
studying a neutral triplet Higgs sector, it is enough to limit ourselves to the 
case in which A3 is zero and A4 = A5. The tree- level results of this triplet 
model are sufficient for demonstrating the phenomenology of the higher 
representations. In this case, the neutral doublet and triplet do not mix and 
the neutral mass eigenstates are 

i^? = /^ = ^{V2x'^ + e\ 

Hl = CHx'' + SHcP'\ Hl = ^{V2e-x''), (10) 

where x° = + iX^')/V2- The masses of the neutral scalars are 
m^o = Sc^Aiu^, m^o' = 3s^A2f^, 

mlro = Xav'^, m\o = 3(A5s|^ + \4Cff)v'^ . (11) 

3 5 

The lightest neutral scalar can be Hf if either ch or Ai is small - the m/y/Xi 
values are shown in Fig. ([1]). It should be noted that charged scalars exist 
with the same masses as and H^. 

The couplings of the neutral scalars to the fermions and the gauge bosons 

are 

^ 2my/CH ' ^ 2mwCH ' ^ 2mwCH ' 

H^W+W- : gmwCH, H^ZZ : mwCH. 

cos"^ 6w 

V<J cos^ OwV'i 
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w- 



V3 



gmwSH, 



H^ZZ 



2ff 



COS' 



mwSH- 



(12) 



It is clear that, at tree-level, the coupling of the to fermions is always en- 
hanced by the factor oH/ch- Conversely, the coupling of the to fermions 
is either enhanced or suppressed, depending on the ratio of sh and ch, and 
the other neutral scalars do not couple to fermions. Importantly, the gauge 
boson couplings to Hf are suppressed by a factor ch with respect to the SM 
and the role of vector boson fusion mechanism for production is reduced 
if Cff is small. The other Higgs couplings to gauge bosons are suppressed by 
Sh- The Hqq couplings presented in Fig. [2] (a) are considerably enhanced 
for small ch in comparison to the SM prediction. In Fig. [2] (b), the HVV 
couplings are shown, again normalized to the SM couplings. 




0.2 0.3 



Figure 2: Couplings of the Higgs bosons to fermions (a) and gauge bosons 
(b), normalized by the Standard Model couplings. 



3.2 Decays of and constraints on the parameters 

In this section, we consider the Higgs boson, which becomes the Standard 
Model Higgs boson for vanishing doublet-triplet mixing. The mass limits for 
Hi can be deduced from the LEP results. The couplings of the Hi to the 
gauge bosons are smaller than in the SM, leading to reduced production of 
Higgs bosons in Higgsstrahlung process [H], through which the Higgs was 
expected to be produced at LEP. The Higgs boson branching ratio to hh 
is 74% for niH = 114 GeV in the Standard Model. Together with other 
fermions, the hh decay mode gives the main contribution to the total width 
of the Higgs boson, and, thus, the Higgs branching ratio does not change 
very much with ch- for mn = 114 GeV the branching ratio changes to 80% 
(81%) for Ch = 0.5 {ch = 0.2). If the Higgs boson is lighter, the change 
is less. If we assume that the number of 6-quark pairs gives the Higgs 
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boson mass limit, it must be heavier than 73 GeV (40 GeV) for ch = 0.5 
{ch = 0.2). Unitarity further constrains most masses, requiring them to be 
less than of the order of 1 TeV [231125]. 

The Yukawa couplings are constrained by perturbativity, which limits 
the Hf coupling to top, 

< V4^. (13) 

2'mwCH 

From this it follows that ch > 0.2, which currently is the most stringent limit 
for cjJl. The latest 95% confidence limit for the cross section x branching 
ratio oi pp ^ H ^ TT for a 120 GeV Higgs boson is observed by the DO 
Collaboration to be approximately 5 pb given 2.2 fb~^ of data [2^. The 
CDF Collaboration observe a similar limit with 1.8 fb~^ of data. Although 
this offers no sensitivity to a SM Higgs boson, which has a production cross 
section of ~1 pb and a branching ratio of ~7%, one would expect increased 
sensitivity for the Hf in the triplet model for small ch (the production 
cross section is increased by a factor of 1/cjj w.r.t to the SM as the main 
production channel would be gluon-gluon fusion via a top quark loop). Thus 
it seems likely that the final combined CDF /DO results will further constrain 
the lower ch limit given the expected integrated luminosity of 8 fb~^ per 
experiment. 

When calculating the branching ratios, it is necessary to consider also 
the loop induced decays of the Higgs bosons to gluons and photons. As 
mentioned, the tree-level gauge boson couplings to Hi are suppressed by a 
factor Ch- The 'yyHf coupling is more complicated. In the SM, the W-loop 
gives the dominant contribution to the jjH coupling. In the triplet model 
however, the fermion coupling is enhanced and the W-coupling suppressed 
for small values of ch, so that the top loop becomes more important and 
the 77-??! coupling is enhanced for small ch- However, since the total width 
increases when ch decreases, the branching ratio to photons remains smaller 
than in the SM. The gluon coupling to is enhanced by 1/ch due to the 
fermion loop. These effects are seen in Fig. [Sj where the branching ratios of 
Hf are presented for rriffO = 120 GeV and 150 GeV. 

4 Central Exclusive DifFractive Production of the 
Triplet Higgs Boson 

The central exclusive production (CEP) of a Higgs boson is defined as pp 
p (B H (B p, where the © denote the presence of large rapidity gaps between 

^Constraints from low energy precision measurements have been obtained from Zbb 
vertex, meson-antimeson mixing and ratios of 6 — » it to 6 — ^ c decays [2^. The radiative 
corrections to the Zbb vertex give the strongest constraint. Assuming rriHs ~ 1 TeV, one 
finds cos 6h > 0.3 with 99.9 % C.L. However, since we consider only tree-level results in 
this work, we have not used this bound. 
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Figure 3: Branching ratios of }i\ to the Standard Model particles for mn 
120 GeV (a) and mu = 150 GeV (b). 



the outgoing protons and the decay products of the central system. It has 
been suggested in recent years that CEP offers a unique complimentary 
measurement to the conventional Higgs search channels, see for example, 
[3 E] , [28] - [33] . Firstly, if the outgoing protons scatter through small angles 
then, to a very good approximation, the primary active di-gluon system 
obeys a = 0, CP-even selection rule |3l]. Here Jz is the projection of the 
total angular momentum along the proton beam axis. The observation of 
the Higgs boson in the CEP channel therefore determines the Higgs quantum 
numbers to be J^^ = 0+"*". Secondly, because the process is exclusive, all 
of the energy /momentum lost by the protons during the interaction goes 
into the production of the central system. Measuring the outgoing proton 
allows the central mass to be measured to just a few GeV, regardless of the 
decay products of the central system. A mass measurement of this type 
will require new forward proton detectors to be installed at ATLAS and/or 
CMS, which we discuss further in section [5Tl 

For a Standard Model Higgs boson, central exclusive diffraction could 
allow the main decay channels (66, WW and rr) to be observed in the same 
production channel, which provides the opportunity to study the Higgs cou- 
pling to 6-quarks. This may be very difficult to access in other search chan- 
nels at the LHC, despite the fact that H ^ bb is by far the dominant decay 
mode for a light SM Higgs boson. Furthermore, CEP can provide valuable 
information on the Higgs sector of MSSM, NMSSM and other popular BSM 
scenarios [35l [36l [311 [33l [STJ [32]. In this paper, we propose that CEP is 
also beneficial if higher representations of the Higgs sector are realized, in 
particular, in searches for the Higgs triplets discussed in section [31 

The theoretical formalism [38] - [l2] for central exclusive production 
contains distinct parts, as illustrated in Fig. [H The cross section can be 
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Figure 4: A symbolic diagram for the central exclusive production of a Higgs 
boson H. 



written in the form [381 [7] 



a{pp ^ p + H + p) 



dQl 

B'' "7 Qi 

(14) 

where B/2 is the t-slope of the proton-Pomeron vertex, (5^) is the soft- 
survival probability and the normalization, N, is given in terms of the H 
gg decay width. The amplitude-squared factor, |...| , can be calculated in 
perturbative QCD because the dominant contribution to the integral comes 
from the region Ag^^^, ^ <C m'jj for the large Higgs mass values of 
interest. The probability amplitudes, fg, to find the appropriate pairs of t- 
channel gluons {xi,x'i) and {x2,X2) are given by skewed unintegrated gluon 
densities at a hard scale /x ~ mH/2. It is important to emphasize, that these 
generalized gluon distributions are usually taken at pt = 0, and then the 
"total" exclusive cross section is calculated by integrating over the transverse 
momentum, px, of the recoil protons. Assuming an exponential behaviour 
results in 



dp^e-'^PT = i/B = (p^). (15) 

Thus, the additional factor in Eq. is not just the gap survival but 
rather the factor {S"^) / B"^ OSS], which has the form S'^ipl)'^ and is much 
less dependent on the parameters of the soft model [121 HIl HI] 

The production cross section for Higgs bosons produced by gluon-gluon 
fusion and decaying to hb is proportional to 

(16) 
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CH = 0.2 


6.35x10"^ 


0.80 
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0.75 


CH = 0.5 


1.01x10^2 


0.79 


1.95x10-3 


0.63 


CH = 0.8 


3.97x10-'' 


0.74 


7.63x10-3 


0.32 


SM 


2.49x10"^ 


0.68 


4.79x10-^ 


0.18 



Table 1: The Hi partial decay width to gluons (expressed in GeV) and 
branching ratio to bb for specific values of Higgs mass and ch ■ The Standard 
Model (SM) prediction is shown for comparative purposes. 



where T{H — > gg) is the decay width to gluons and BR(ff — > 56) is the 
branching ratio to bb quarks. Table [1] shows the value of these parameters for 
the SM Higgs and the lightest Higgs, Hf, in the triplet model. The central 
exclusive H ^ bb cross section can therefore be enhanced by a large factor 
with respect to the Standard Model - we discuss this further in Section [5] 

The CEP formalism has been extensively checked using the difFractive 
production of J/ip and the leading neutron spectra [l5] at HERA and the 
CDF data on central exclusive production processes [l6] - [IH]. Further 
tests of the formalism using the early LHC data have also been suggested 
[42| . The main uncertainties are associated with: 

• The probability (S^) that additional secondaries will not populate the 
gaps. 

• The probability to find the appropriate gluons that are given by gen- 
eralized, unintegrated distributions fg{x,x' ,Qf). 

• Higher order QCD corrections to the hard subprocess, in particular, 
the Sudakov suppression. 

• The so-called 'enhanced absorptive corrections' [HI [501 131] a-nd other 
effects that may violate the soft-hard factorization. 

Let us focus first on the gap survival factor (5"^). Since soft physics 
is involved, we need a reliable model of soft interactions to quantify the 
role of the absorption effects. In |39t HUl [JT] soft interaction models were 
developed and tuned to describe all the available high energy soft pp in- 
teraction data. These models account for (i) elastic rescattering (with the 
two protons in intermediate states), (ii) the probability of low-mass pro- 
ton excitations, and (iii) the screening corrections due to high-mass proton 
dissociation (enhanced absorption). The first two effects result in the so- 
called eikonal contribution, (5'gj]^). In the most recent version of the soft 
rescattering model [111133], the KMR group obtained (5'^j;^)cfT = 0.025 when 
adjusting (S^jj^^) to its value corresponding to an exponential slope B = A 
GeV-2. 
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In the presence of enhanced screening, however, there is no longer exact 
factorisation between the hard and soft parts of the process, see for example 
[50l HSl H]. The latest calculations [IT] indicate that, in the case of the 
SM Higgs production at the LHC, the effective survival factor, due to both 
eikonal and enhanced rescattering, is (S'^)ofr — 0.015 it 0.01. It should be 
noted that the exclusive dijet, 77 and Xc production data from CDF and 
the leading neutron data at HERA indicate that {S^^^^ is somewhat larger 
than this estimate, such that {S'^)es is nearer the upper limit of the quoted 
interval. In any case, it will be possible to measure {S"^^^) using early LHC 
data 021 EH. 

As the generalized, unintegrated gluon distribution fg has not been mea- 
sured explicitly, it is obtained in the KMR approach |38l |7j from the conven- 
tional gluon distribution, g{x,Q^), known from the global parton analyses. 
The main uncertainty comes from the lack of knowledge of the integrated 
gluon distribution at low x and small scales. It was found in [12] that a 
variety of recent global analyses give a spread of 

xg = (3-3.8) for x = and xg = (3.4-4.5) for x = 10~^(17) 

for = 4 GeV^. These are big uncertainties bearing in mind that the 
CEP cross section depends on (xg)^. A similar estimate of the uncertainty 
from the input gluon distribution functions was presented in [33j. The un- 
certainty related to the Sudakov factor is addressed in [32], along with the 
measurements that could reduce the uncertainty using early LHC data. 

The overall uncertainty factor in the calculation of the CEP of Higgs 
bosons at the LHC was estimated to be approximately 2.5 in [351131] . Again, 
we note that the first LHC runs will allow the accuracy to be drastically 
improved. 

5 Simulation of Higgs production in the triplet 
model 

5.1 Forward proton tagging 

The forward proton detectors will need to be installed in the high disper- 
sion regions 220 m and 420 m either side of the interaction point at ATLAS 
and/or CM^. We restrict our discussion and analysis to the ATLAS in- 
teraction point (IP), and note that a similar result would be obtained at 
CMS. There are three important aspects of forward proton tagging at the 
LHC that need to be considered for the purposes of this analysis; the ac- 
ceptance and resolution of the proposed forward proton detectors and the 
ability of the detectors to measure the time-of-flight of each proton from the 
interaction point. 

®We refer the reader to the FP420 R&D report for more details on the detectors [8]. 
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The acceptance of the forward proton system depends on the distance 
that each active detector is from the LHC beam. We choose the proton 
detectors located at 220 m (420 m) from the IP to be 2 mm (5 mm) from 
the beam and use the FPtrack program to track the path of the protons 
through the LHC lattice in order to fully determine the acceptance. The 
acceptance of the forward proton detector system is dependent on the mass 
of the centrally produced object, M, which is given by 

M2 = eie2s (18) 

where and .^2 are the fractional longitudinal momentum losses of the 
outgoing protons. For central masses less than 200 GeV, the protons are 
either both detected at 420m (symmetric tagging) or one is detected at 220 
m and one at 420 m (asymmetric tagging). For example, the acceptance 
for a central mass of 120 GeV is approximately 28% for symmetric events 
and 20% for asymmetric event^ [8j. This acceptance changes to 20% and 
40% respectively for a 150 GeV central mass. Furthermore, symmetric and 
asymmetric events also have different mass resolution for a given central 
mass; the resolution of a 150 GeV Higgs boson is approximately 1.5 GeV 
and 4 GeV for symmetric and asymmetric proton tagging respectively. 

Finally, the forward proton detectors will be capable of measuring the 
time-of-flight of each proton from the interaction point to an accuracy of 
10 ps. The difference in the time-of-flight measurement of the protons gives 
a measurement of the interaction vertex to 2.1 mm, assuming that the ref- 
erence timing system has negligible jitter. This vertex reconstruction proves 
to be very useful in background rejection, as discussed in Section [531 

5.2 Signal and background event generation 

The central exclusive signal and background events are simulated with full 
parton showering and hadronization effects using the ExHuME vl.3.4 event 
generator |53j . ExHuME contains a direct implementation of the KMR 
calculation [38^ \7\ of central exclusive diffraction given in Sec. HI The 
CTEQ6M [51] parton distribution functions are used to calculate the gener- 
alized gluon distributions, fg. The generator level cross sections for central 
exclusive H ^ bb production in the triplet model are presented in Table [2] 
for niH = 120, 150 GeV and ch = 0.2, 0.5, 0.8. In Section [531 we will present 
results for each of these parameter choices. 

The backgrounds to H ^ bb can be broken down into three broad cate- 
gories; central exclusive, double pomeron exchange and overlap. We use the 
ExHuME event generator for the central exclusive backgrounds. ExHuME 
contains the leading order calculation for central exclusive bb production. 

^Dead material in the detectors at 220 m can reduce the acceptance for symmetric 
tagging. We do not consider that effect here. 
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m// = 150 GeV 


CH = 0.2 


113.5 


55.2 


CH = 0.5 


18.0 


7.4 


CH = 0.8 


6.6 


1.5 



Table 2: Generator level cross sections, CTj^^f^^, for central exclusive Higgs 
boson production for mn = 120, 150 GeV and ch = 0.2, 0.5, 0.8. 

Recent results however |55j . show that the one- loop corrections to this pro- 
cess reduce the cross section by a factor of approximately two. We normalize 
the bb events generated by ExHuME accordingly. Central exclusive gg pro- 
duction has a much larger cross section than bb and can act as a background 
when the gluon jets are mis-identified by the b-tagging algorithms. At AT- 
LAS, the mis-identification rate for each gluon jet is 1.3% for a b-tagging 
efficiency of 60% |56| . We do not generate exclusive cc events for two reasons. 
Firstly, due to the Jz = spin-selection rule ^34J, exclusive cc production 
is suppressed with respect to bb by a factor of m'^/ml. Furthermore, the 
mis-tag rate for c-jets to be identified as b-jets is ~10%. Thus, the back- 
ground contribution from exclusive cc events is considered to be negligible 
in comparison to the exclusive bb events. Higher order events, such as bbg 
have been studied in [57] with the conclusion that these types of events 
should be negligible after all experimental cuts. The demonstration of this 
using event generators cannot be completed until the relevant processes are 
implemented into the ExHuME Monte Carlo. 

Double pomeron exchange (DPE) is the process pp p+X+p, where the 
central system, X, is produced by pomeron-pomeron fusion. The pomeron 
is assigned a partonic structure and so there are always 'pomeron remnants' 
accompanying the hard scatter. DPE has been extensively studied in rela- 
tion to H ^ bb and it has been concluded that this background is negligible 
after appropriate experimental cuts [331 [58]. We do not consider these types 
of events further. 

In addition to these standard backgrounds, we also examine the effect of 
the overlap backgrounds. This source of background is important when there 
are a large number of pp interactions in each bunch crossing at the LHC. 
The largest overlap background is a three-fold coincidence between two soft 
events {pp — > pX), which produce forward protons within the acceptance 
of the forward detectors, and an inelastic event, which produces the hard 
scatter pp ^ X and, thus, can mimic our signal. To simulate these events 
we use HERWIG [59] plus JIMMY [60] to generate pp bb, using the tune 
(A) to Tevatron data [HI]. The forward protons are then added into the 
event using the prescription given in [33], which also allows us to calculate 
the probability of the coincidence as a function of the LHC luminosity. 

The overlap background is initially reduced using the vertex matching 
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provided by the proton time-of-flight (TOF) information. As the protons 
do not come from the same interaction as the jets, the vertex reconstructed 
using TOF will not, in general, coincide with the di-jet vertex. Given a fast- 
timing resolution of 10 ps, a rejection factor of 18 (15) can be obtained at 
low (high) luminosity by requiring that the di-jet vertex be within ±4. 2mm 
of the 'fake' vertex reconstructed from TOF. Approximately 95% of the CEP 
events will be retained by this requirement. 

We do not consider the backgrounds from two-fold coincidences. It was 
demonstrated in [33] that the largest of these backgrounds - the coincidence 
between a soft central diffractive scattering and a standard QCD 2 — > 2 scat- 
tering - was at least a factor of five smaller than the threefold coincidence. 
Furthermore, as discussed in the FP420 R&D report [8], this background 
is (i) probably overestimated and (ii) will be additionally rejected by the 
charged track cut outlined in the next section. 

When generating the event samples, we require that the central mass 
be in the range 80 < M < 250 GeV, which improves the event generator 
efficiency and is the broad region of interest for this study. To approximate 
the detector effects, we smear the energy, momenta and angles of each central 
final state particle according to the ATLAS detector resolution [62]. The 
outgoing forward proton momenta are smeared by the amount specified in 
[63] . A mid-point cone algorithm is then applied to the samples and events 
retained if the leading jet has transverse energy greater than 45 GeV. Finally, 
b-tagging efficiencies are imposed after matching the two leading jets to the 
partonic level. 

5.3 Experimental cuts 

To enhance the signal, we follow the experimental method used in a previous 
study of H ^ bb in the SM and the MSSM [3H], which imposes a number of 
exclusivity cut^. Firstly, the rapidity of the central system can be estimated 
from the forward proton detectors by 



The difference. Ay, between this rapidity measurement and the average 
pseudo-rapidity of the di-jets should be approximately zero for an exclusive 
event, Exclusive candidates satisfy 



where rji and 772 are the pseudo-rapidities of the leading jets. 

*A somewhat different experimental method was discussed in 3f . ffowever, the final 
experimental efficiency broadly agrees with that used here [8^. 




(19) 




(20) 
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The di-jet mass fraction, Rj, determines the fraction of the central mass 
that is contained within the di-jet system. Rj is an improved version of 
the Rjj variable [64J, which was used by CDF in the search for exclusive 
di-jet events [46j. For an exclusive di-jet event, one expects all of the mass 
to be contained in the di-jets, and hence Rj = 1. However, parton show- 
ering/hadronisation effects can result in energy outside of the jets. Fur- 
thermore, detector resolution will smear the di-jet mass measurement. An 
exclusive event is defined to be 

0.85 < i?, = ^cosh(ryi-y) < 1.1, (21) 

where Ej^ is the transverse energy of the leading jet. 

The third exclusivity cut requires the di-jets to be back to back in az- 
imuth, which reflects the suppression of initial state radiation for an exclusive 
di-jet event. The back to back requirement is 

Ivr - A(/)| < 0.15, (22) 

where A(p is the azimuthal angle between the jets. It is also possible to 
examine each event for underlying event activity, caused by multiple parton- 
parton interactions. The exclusive events do not have these additional scat- 
ters because the proton remains intact. It is possible to reject inclusive 
events, and, hence, the overlap background, by requiring few charged tracks, 
Nc, associated with the di-jet vertex but outside of the jets. This definition 
is, of course, dependent on the jet algorithm used to define the jets. An 
algorithm independent approach is to examine the charged track activity 
perpendicular to the leading jet, A*"^. In this approach, charged tracks are 
assigned to the underlying event if they satisfy 

TT , , , , 27r 47r , , , , 57r , , 

3 < |'/'fc-</'i| < Y or —<\cPk-cPi\<—, (23) 

where (f)^ is a the azimuthal angle of a charged track and (f>i is the azimuthal 
angle of the leading jet. In this analysis, we identify exclusive events by 

A^c < 3 and < 1. (24) 

For completeness, in Table [3] we present the final cross sections for the 
signal {niH =120 GeV, ch = 0.5) and background events. It should be noted 
that the signal is concentrated at M = 120 GeV, whereas the backgrounds 
form a continuum across the mass range 80<M<250 GeV. For details on 
the efficiencies of the individual cuts, and motivation for the cut choices, we 
refer the reader to [33]. The overlap backgrounds are luminosity dependent 
and are presented for constant luminosities of 10^^ cm~^ (low) and 
10^^ cm~^ s"-*^ (high). 



15 



Generator 


Process 


^■420-420 (fb) 


^■420-220 (fb) 


ExHuME 


H ^bb 


0.53 


0.28 




bb 


0.53 


0.27 




99 


1.08 


0.91 


Overlap (L) 


bb 


0.07 


0.09 


Overlap (H) 


bb 


11.0 


13.7 


Total bgrd (L) 




1.68 


1.27 


Total bgrd (H) 




12.6 


14.9 



Table 3: The final cross sections for the H ^ bb {mn = 120 GeV, ch = 0.5) 
and relevant backgrounds after the all cuts discussed in the text. The overlap 
backgrounds are defined at both low (L) and high (H) luminosity. All the 
backgrounds form a continuum over the range 80<M<250 GeV. 

5.4 Trigger strategies 

A major experimental challenge for central exclusive jet analyses is devel- 
oping a trigger strategy to retain enough events. At ATLAS, jets with 
Et ~ 50 GeV are heavily prescaled in the level one (LI) trigger, due to the 
high rate and the lack of additional rejective power in the high level trigger 
(HLT). The total LI rate allowed at ATLAS is 75-100 kHz, which must be 
reduced to ~100 Hz after the HLT. In this section, we discuss three possible 
trigger strategies. The first possibility is to exploit the muon rich nature 
of bb events. The lowest muon threshold (MU6) at ATLAS is designed to 
retain 80% of muons with transverse momentum greater than 6 GeV. The 
single muon trigger efficiency for bb events is 11% [33]. In order to keep the 
LI rate down it will be necessary to require that the event contains jet with 
Et > 40 GeV in conjunction with the muon. 

The second trigger strategy is to require a 40 GeV jet in conjunction with 
a proton tagged in a detector at 220 m from the interaction poinl|£|. This 
trigger has been extensively studied in previous work and it is expected 
that the unprescaled LI rate will be less than 1 kHz up to a luminosity of 
L=2xl0^^ cm^^ s~^. This rate however, will scale with L^ and the trigger 
may have to be prescaled to give a fixed rate at the highest luminosity. We 
will investigate two fixed rate triggers; R5 is a LI rate of 5 kHz and RIO 
is a LI rate of 10 kHz. The drawback of this trigger strategy is that the 
symmetric events will not be retained, which could potentially have a large 
impact for a light Higgs boson measurement. 

It will be possible to dramatically reduce any high LI rate using ad- 
ditional information in the HLT. Firstly, requiring that there be a proton 
detected at 420 m and using time-of-flight information to match the vertices 

®The information from the detectors at 420m will not reach the central trigger processor 
within the latency of 2.5/is and so cannot be used in the LI decision. 
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will reject the events by a factor of approximately 600 (60) at low (high) 
luminosity. A loose b-tagging requirement could also be added. Finally, the 
Rj and Ay cuts reject the overlap events by a factor of approximately 100 
and could be used to further reduce the rate. 

The final trigger strategy is to allow a larger LI rate for the 40 GeV jets, 
which is then reduced in the HLT by requiring two in-time protons. This 
trigger has also been studied in previous work [33]. The advantage of this 
approach is that it would retain both symmetric and asymmetric events. 
The LI rate for 40 GeV jets is expected to be approximately 25 kHz at 
low luminosity, rising to 250kHz at high luminosity. We define a trigger for 
40 GeV jets, JR25, which has a fixed rate of 25 kH^^. This trigger would 
be unprescaled at low luminosity and prescaled by a factor of ten at high 
luminosity. 

5.5 Significance of observation and expected mass distribu- 
tions 

In this section, we estimate the significance of observing a neutral Higgs 
boson in the triplet model for the parameter choices presented in section [5^2] 
and for the trigger strategies outlined in section 15.41 As the overlap back- 
ground is luminosity dependent (Table [3]), we must specify how the data was 
collected. For example, we examine the significance for an integrated lumi- 
nosity of 60 fb~^, which corresponds to between three and four years of data 
acquisition given a peak luminosity of 2x10'^^ cm~^ s~^. We also present 
results for 300 fb~^ of data, which corresponds to between three and four 
years of data acquisition given a peak luminosity of 10^^ cm~^ s~^. It should 
be noted, however, that the data acquisition at the LHC will not be collected 
at one specific luminosity. Firstly, the peak luminosity will increase during 
the lifetime of the LHC with an improved understanding of the machine. 
Secondly, during each store, the luminosity will decrease exponentially with 
a lifetime of approximately 14 hours. We crudely approximate these effects 
by defining that half the data is collected at the peak luminosity and half 
the data at 75% of the maximum. 

The signal and background events that pass the selection criteria are 
normalized to an expected number, A^, for each trigger /luminosity scenario, 
i.e. 

N = Lae, 

where L is the total integrated luminosity, a is the final cross section for 
each process as shown in Table [3] and e is the trigger efficiency. In the case 
of the overlap background, a is dependent on the assumed peak luminosity 

^'^It was shown in CMS-based study [66] that this rate can be reduced by a factor of 
two by requiring that the majority of the transverse energy in the detector be contained 
within the dijets - i.e. that {E^ + E^) / Ht > 0.9, where Ht is the scalar sum of transverse 
energy deposited in the detector [65j . 
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as described above. The MC distributions are then used to predict the 
expected number of events in each bin of a mass distribution. A pseudo- 
experiment is then carried out by randomly picking a number of events 
for each bin in the mass distribution according to a Poisson distribution. 
The significance of each pseudo-experiment is then obtained by fitting the 
pseudo-data with a signal-plus-background function and a background-only 
function. The significance is estimated from the difference in of the two 
fits by 



We assume the background function will be well known from data - in this 
analysis we use the weighted MC background events to determine the shape 
and allow the normalization to var}0. The signal function is a Gaussian and 
all parameters are allowed to vary. We repeat the procedure for 500 pseudo- 
experiments to determine the average significance of each luminosity/trigger 
scenario and ensure that our results are consistent and the presented distri- 
butions typical. 

To determine the significance for each parameter choice, the trigger 
strategies outlined in Sec 15.41 are evaluated and the best method chosen 
to retain the events. We choose the mn = 120 GeV, ch = 0.5 point as 
our reference. Figure [5] (a) shows the significance as a function of luminos- 
ity for the R5, RIO and JR25 trigger strategies given three years of data 
acquisition at that luminosity. It is clear that the JR25 trigger is the best 
choice at low luminosity as it retains a high fraction of both symmetric 
and asymmetric events. At higher luminosities, the R5/R10 triggers be- 
come more favourable; the R5 (RIO) trigger do not become prescaled until 
4.5x10^^ cm~^ (6.3x10'^^ cm~^ s~^) and the number of exclusive events 
in the final sample is increased even though the symmetric events are not 
retained. Figure [5] (b) shows the significance for R5, RIO and JR25 when 
combined with the MU6 trigger. 

Figure [6] shows the expected mass distributions for niH = 120, 150 GeV 
and Ch = 0.2,0.5 for 60 fb~^ of data collected at 2x10'^^ cm^^ s~^ using 
a JR25+MU6 trigger. Example distributions for ch = 0.8 are not shown 
as the calculated significance is less than 3a. Figure [7] shows the same 
distributions but for 300 fb~^ of data collected at lO^'' cm~^ s~^ using 
the R10+MU6 trigger. A summary of the significance for each parameter 
choice is presented in Table H] for both luminosity/trigger scenarios. The 
significance is approximately 4o" for the ch = 0.5 parameter points and 12a 
for the Ch = 0.2 parameter points. 

For each parameter point, the information obtained from the fitting pro- 
cedure in each pseudo-experiment can be used to obtain an RMS spread of 
Higgs masses. This gives a reasonable estimate of the error on the Higgs mass 

^^We have checked our resuhs using a quadratic background function and observe httle 
difference. 
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Figure 5: (a) Expected significance {mn = 120 GeV, ch = 0.5) for the R5, 
RIO and JR25 triggers as a function of luminosity given three years of data 
acquisition at each luminosity. We assume that 10 fb~^ of data is collected 
each year per 10^^ cm~^ of luminosity, i.e. 120 fb~^ of data is collected 
at 4x10^3 cm s . (b) The significance of the R5, RIO and JR25 triggers 
when combined with the MUG trigger. 



as measured by the forward proton detectors. Table [S] shows the RMS spread 
for each (significant) parameter choice. The RIO trigger strategy (used at 
high luminosity) retains only events with asymmetrically tagged protons, 
whereas the JR25 trigger (used at low luminosity) retains both symmetric 
and asymmetric events. Thus the high luminosity scenario in Table [5] results 
in a poorer mass measurement than the low luminosity scenario because of 
the lower fraction of symmetric events in the sample. The final trigger strat- 
egy at ATLAS/CMS will have to balance the need for observation with the 
opportunity for more precise measurements. For all parameter choices the 
mass measurement can be made to better than 2 GeV; for ch = 0.2 the 
mass measurement is better than 0.3 GeV. 

It should be noted that the significances obtained for the ch = 0.2 
parameter points are well in excess of lOcr and the high rate triggers assumed 
in this analysis are not strictly needed for the measurement. Indeed, the 
MUG trigger alone is capable of retaining enough events for the analysis 
- the significance for a 120 GeV Higgs boson is 4.5 for GO fb~^ collected 
at 2x10'^^ cm^^ s^^. The disadvantage is that the mass measurement is 
somewhat degraded, due to the reduced number of events. 

The Ch = 0.8 points are observable if the overlap background can be 
additionally rejected. There are two possibilities related to improvements in 
the fast-timing system. Firstly, if the time-of-flight resolution of the forward 
detectors is improved, then both the vertex resolution and overlap rejection 
is improved by the same factor. Secondly, if the central jets can be timed 
to an accuracy of 70 ps, using optimal signal filtering in the ATLAS Liquid 
Argon Calorimeter [57], a factor of two rejection on the overlap background 
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(c) (d) 

Figure 6: Expected mass distributions given 60 fb^^ of data, collected at 
2x10'^^ cm~^ s~^ using a JR25+MU6 trigger, for the following parameter 
choices: (a) uih = 120 GeV and ch = 0.2, significance is 12.7a (b) mn = 
150 GeV and ch = 0.2, significance is ll.Do". (c) tuh = 120 GeV and 
Ch = 0.5, significance is 4.5(j. (d) mn = 150 GeV and ch = 0.5, significance 
is 3.9a. 



would be gained. If the central timing could be performed at the 10 ps 
level, then a factor of 12 rejection would be observed. The significance of 
the Ch = 0.8 parameter point for a 120 GeV Higgs boson increases to 3.2a" 
at high luminosity if the overlap background is additionally rejected by a 
factor of five. Improved overlap rejection would also increase the significance 
of the other parameter points at high luminosity. 



6 Doubly charged Higgs bosons 

An important feature of the triplet model is the existence of doubly charged 
Higgs bosons. Although not the focus of this paper, we comment on the 
possibility to observe these using the forward proton system discussed in 
section 15.11 The process of interest is — > p © H^^H © p, where the 
central system is produced via photon-photon fusion, i.e. 77 H++H—. 
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Figure 7: Expected mass distributions given 300 fb ^ of data, collected at 
10^4 cm ^ s ^ using a R10+MU6 trigger for the following parameter choices 
(a) ran = 120 GeV and ch = 0.2, significance is 13.7a (b) uih = 150 GeV 
and Ch = 0.2, significance is 12.7cj. (c) ttih = 120 GeV and ch = 0.5, 
significance is 4.3cr. (d) niH = 150 GeV and ch = 0.5, significance is A.3a. 
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ruH (GeV) 


CH 


Significance (a) 






60 fh-^ 


300 fh-^ 


120 


0.2 


12.7 


13.7 




0.5 


4.5 


4.3 




0.8 


2.6 


2.5 


150 


0.2 


11.9 


12.7 




0.5 


3.9 


4.3 




0.8 


2.0 


2.1 



Table 4: Significance of central exclusive Higgs boson measurement for 
niH = 120, 150 GeV and ch = 0.2, 0.5, 0.8. The significance is obtained 
after using the cuts, trigger strategies and fitting procedure outlined in the 
text. As the overlap background is luminosity dependent, we present results 
for both 60 fb~^ of data taken at a peak luminosity of 2x10^^ cm~^ s~^ and 
for 300 fb~^ of data taken at a peak luminosity of 10^^ cm~^ s~^. 



ruH (GeV) 


CH 




(GeV) 






60 fb-i 


300 fb-i 


120 


0.2 


0.2 


0.3 




0.5 


0.8 


4.5 


150 


0.2 


0.2 


0.3 




0.5 


1.8 


2.4 



Table 5: The error in the Higgs boson mass measurement, 0"^^, esti- 
mated from the RMS spread of mass measurements made in the pseudo- 
experiments, for ruH = 120,150 GeV and ch = 0.2,0.5. The 60 fb~^ of 
data (taken at a peak luminosity of 2x10'^^ cm"^ s""*^) uses a JR25 + MU6 
trigger strategy whereas the 300 fb^^ of data (taken at a peak luminosity of 
10^^ cm-2 s"^) uses RIO + MU6. 

As two Higgs bosons are produced, the central system will be of higher mass 
than discussed in previous sections and at least one of the protons will be 
tagged by detectors placed at 220m from the IP. 

The production cross section could be large [68J. It was shown in [691 [8] 
that the production cross section was 0.07 fb for a 200 GeV for singly charged 
scalar pair production. However, the cross section for doubly charged Higgs 
boson pair production is a factor of 16 larger due to the factor of two 
increase in charge. Thus the cross section for pair-production of a 200 GeV 
doubly charged Higgs boson increases to 1.1 fb. It should be noted that the 
photon fusion production cross section decreases roughly by ~ [7], 
and so the pair production of a 300 GeV doubly charged scalar would be 

^■^The doubly charged Higgs boson mass in the Machacek-Georgi model is the same as 
the -fff scalar mass, and therefore, JJ""*"^ with mass m^++ ~ 200 — 300 GeV is possible, 
as it is heavier than the lightest neutral Higgs H^. 
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approximately 0.3 fb. 

The forward proton acceptance for central systems in the range 400 GeV- 
1 TeV is better than 40% [63j . This means that we would expect more than 
130 events with both outgoing protons measured for the pair production 
of a 200 GeV doubly charged Higgs boson given 300 fb^^ of data acquisi- 
tion. This decreases to 40 events if the mass of the doubly charged Higgs is 
300 GeV. 

The events should be retained using the standard electron/muon trigger 
strategies currently in place at ATLAS/CMS. A doubly charged Higgs that 
is heavier than 200 GeV has several possible decay modes: H~^^ W~^W~^, 
W~^H^, H~^H~^ , . In any case, a high fraction of events will contain at 
least one electron/muon in the final state. 

Therefore, forward proton tagging allows the possibility to study doubly 
charged scalars that are lighter than around 300 GeV. Note that, in exclusive 
production, the background conditions are more favourable in comparison 
to the conventional pp H^^H X case considered e.g. in [68] . Also, 
in principle, the forward proton mode may allow to measure the mass 
more accurately than in the inclusive case. 

7 Conclusions 

Searches for the manifestation of the extended Higgs sector at the LHC may 
allow new insight in the nature of electroweak symmetry breaking. The cen- 
tral exclusive production mechanism would provide a very powerful tool to 
complement the standard strategies at the LHC for studying Higgs particles. 
Here we focus on the production of the neutral Higgs boson of the triplet 
model in the forward proton mode. We assume a model with the tree-level 
value of the electroweak p-parameter consistent with experiment, p = 1. 
Although this model is used as a benchmark model for the triplets, our re- 
sults are more general. An extra contribution from other representations 
enhances the doublet Yukawa couplings resulting in a different experimental 
signature to that of the SM. We show that a factor of two enhancement 
of the fermion couplings due to the higher representations implies a signif- 
icant difference to the doublet case. Let us emphasize that in the case of 
the current model, all the fermion couplings to the Higgs boson, which is 
responsible for the fermion masses, increase. This is in contrast with, for 
instance, the MSSM, where couplings of up-type and down-type fermions 
change from the Standard Model differently, due to the fact that there are 
only doublets in the model. It is a common feature of higher Higgs represen- 
tations that the doublet couplings are enhanced, which thus indicates that 
higher representations are involved. 

We present a detailed Monte Carlo analysis of the central exclusive pro- 
duction of a triplet model Higgs boson for a number of parameter choices. 
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For ch < 0.5, we have shown that a hght Higgs boson (of mass 120- 
150 GeV) can be observed with a 4c7 (or better) significance if a fixed rate 
trigger is used. We find that a fixed rate single jet trigger is optimal at 
low luminosities whereas a fixed rate forward proton trigger (i.e. one pro- 
ton detected at 220m in conjunction with a central jet) is optimal at high 
luminosities. 

The expected error in the Higgs mass measurement using forward proton 
detectors is small. For ch = 0.2, we find that the mass of the Higgs boson is 
measured to better than 0.3 GeV, regardless of the luminosity/trigger sce- 
narios. This is due to the excellent mass resolution of the forward detectors 
and the large number of events. For ch = 0.5, there are less events and 
so the error in the mass measurement increases. However, regardless of the 
parameter choice, the mass measurement can always be made to better than 
2 GeV if a fixed rate single jet trigger is used to retain events in which both 
protons are tagged at 420 m from the IP. 
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